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Abstract 

This paper describes a model for magnetization reversal in polycrystalline Ferro- 
magnetic/ Antiferromagnetic exchange biased bilayers. We assume that the exchange 
energy can be expanded into cosine power series. We show that it is possible to fit 
experimental asymmetric shape of hysteresis loops in exchange biased bilayer for 
any direction of the applied field. The hysteresis asymmetry is discussed in terms of 
energy considerations. An angle [3 is introduced to quantify the easy axis dispersion 
of AF grains. 
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1 Introduction 



Exchange bias is a phenomenon known for a long time [1] and characterized 
by a displacement (of shift value ifg), often called "bias field", of the hysteresis 
loop along the field axis as a result of exchange coupling between ferromag- 
netic (F) and antiferromagnetic (AF) layers. The original model proposed by 
Meiklejohn and Bean [1] is based on a perfect uncompensated monolayer of 
spin at the surface of the antiferromagnetic layer. This model predicts a shift 
field Hg that is two orders of magnitude larger than that observed experi- 
mentally. Several recent theories give improved predictions of the magnitude 
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of Hs, but they interpret in a different way the physical origin of the effect. 
For example, Mauri et al. [2] have shown that realistic values for Hg can 
be obtained if one assumes a domain wall formation in the AF layer during 
the reversal of the F magnetization. Using locahzed atomic spins, Koon [3] 
suggested 90 degrees spin-flop couphng between the AF and F layers and pre- 
dicted a correct magnitude oi Hg. A different approach was undertaken by 
Malozemoff [4] who assumed that the interface roughness greatly reduces the 
number of uncompensated spins at the AF surface giving rise to a smaller bias 
field. Siles and McMichael have proposed a model for polycrystalline F/AF 
bilayers [5] consisting of a F layer interacting with independent AF grains 
using three contributions for the energy of each grain, a direct coupling, a 
spin-flop coupling and a partial domain wall in the AF grains. They showed 
that the unidirectional anisotropy comes from stable grains as the magneti- 
zation is rotated while the hysteretic effects come from grains where the AF 
order is switchable. In all cases the spin-flop coupling does not contribute to 
the unidirectional anisotropy. Despite many theoretical studies of the origin 
and the magnitude of the bias field, a complete description of the exchange 
anisotropy should incorporate explanations for the other features of the ex- 
change biased systems such as the directional properties of the exchange bias 
and the strong asymmetry of the hysteresis loop that occurs in these systems. 
The angular dependence of Hg and the coercivity He was first explored experi- 
mentally in NiFe/CoO bilayers [6]. The variation of Hg and He versus the field 
angle 9h was described with cosine Fourier series, with odd and even terms 
for Hg and H^ respectively instead of simple sinusoidal functions as suggested 
initially [1]. The angular dependence of the exchange anisotropy in F/AF 
bilayers have been studied both theoretically and experimentally by several 
authors. For some of them [6,7], the complicated angular behavior comes from 
the larger contribution of higher-order terms a„cos(n 9h) (n > 3) in Hg and 
He- Kim et al. [8] have modeled the F/AF bilayer as a cubic structure and 
used an Heinsenberg Hamiltonien to examine the angular dependence of per- 
fect interfaces. For compensated interface, their numerical results of angular 
dependence of the exchange field were fit following the work of Ambrose et al. 
[6], who have shown that the angular dependence has a simple cosine function. 
For rough interface the behavior is more complicated and simple cosine form 
is suppressed and high order terms appear. Kim et a/. [8] suggested that the 
possible origin of the higher order terms may come from the nonrigid spin 
structure of the AF layer. Geshev et al. [9] assumed a planar domain wall 
formation at the antiferromagnetic side of the interface with the reversal of 
ferromagnetic orientation [9]. Xi and White[10] have shown that the complex 
angular dependence of the exchange bias in NiFe/CrMnPt bilayers can be un- 
derstood by a simple Stoner-Wohlfath model involving only a uniaxial and 
a unidirectional exchange coupling. Very recently Krivorotov et a/. [11] have 
shown that the exchange anisotropy, induced by the F/AF exchange coupling, 
consists of a biaxial component, which gives rise to the enhanced coercivity. In 
the latter model, the unidirectional component is responsible for the shift of 
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the hysteresis loop and the threefold symmetry component for the symmetry 
breaking in the magnetization reversal process. Tang et al. [12] have recently 
shown that the exchange bias effect, described by a unidirectional anisotropy, 
is also accompanied by induced uniaxial and fourfold in-plane contribution in 
Fe/MnPd exchange biased bilayers. Most of these authors obtain good agree- 
ment between theoretical and experimental angular variation of Hg and He- 

Many authors have only fit the angular dependence [6,8]. Xi et al. [10] have 
used the fit parameters for the angular dependence of the coercivity and the 
exchange field to calculate the hysteresis loops for different directions of the 
applied field. For the easy axis hysteresis loop, the coercivity is about two time 
larger than the measured one and the shifted hysteresis loop is symmetric. Our 
approach is different, first we fit the easy axis hysteresis loop and then, by using 
the fit parameters, calculate the angular dependence of the coercivity and the 
shift field and also the hysteresis loop for any angle of the apphed field. 

As to the asymmetry of the hysteresis loops in exchange biased bilayers, there 
is a large variety of asymmetric shapes. The modeling of hysteresis loops 
and especially its asymmetry was achieved only for some particular shapes 
[2,13,14]. However none of these above methods describe satisfactorily the 
shape of the measured easy axis hysteresis loop presented on Fig. 2 (a). This is 
why we propose a theoretical model based on the development of the energy 
functional for the exchange field in cosine powers series. We show that the 
experimental asymmetric hysteresis loop measured along the easy axis can be 
fit with our theoretical model. In the present work, we study the magnetiza- 
tion reversal process using conventional (VSM) magnetometry in NiFe/MnNi 
bilayers. The theoretical magnetization curves, calculated for different orien- 
tation of the applied field, are in good agreement with the experimental data. 
We show that the asymmetry in the magnetization reversal may be understood 
in terms of energetic considerations. 



2 Sample prepeiration and analysis 

Substrate\Ari8iFei9(520A)\Mn^Ariioo-x(800A) bilayers were grown on Corn- 
ing Glass substrate by RF diode sputtering using a standard Z 550 Leybold 
equipment with a magnetic field of 24 kA/m Oe applied during deposition 
to induce an uniaxial anisotropy. The background pressure was lower than 
4 X 10~^ mbar. Ni chips were homogeneously added to a four inches diameter 
Mn target in order to get films in the Mn composition range 5-80 percent. The 
chemical homogeneity was verified by Electron Probe Micro Analysis (EPMA) 
on several points of the sample. The Mn composition variation is about one 
percent on the entire sample. As the as-deposited samples did not exhibit 
exchange bias, after deposition, they were annealed in a magnetic field of 80 
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kA/m, aligned with the easy axis of the film, at 300°C for 5 hours to induce the 
exchange field. The crystallographic structure was examined by X-ray difi'rac- 
tometry with Cu Ka radiation. A (111) texture for the NigiFeiQ was favored. 
Annealed MnxNiiQQ_x films deposited on the underlying NigiFcig film were 
found to have a (111) texture with a fct structure in the composition range 40 
<x<80 that exhibits exchange bias [15]. The magnetic properties, such as the 
saturation magnetization Ms, were obtained from magnetization loops (M-H 
loops) measured at room temperature using a VSM. It is well established that 
the shift of the hysteresis loop of exchange biased bilayers depend on intrinsic 
properties such as the AF layer thickness, the F layer thickness and the depo- 
sition conditions (for a review see [16]). All these parameters induce hysteresis 
loop with a more or less pronounced asymmetry. For our study we have cho- 
sen samples with a well-defined magnetization reversal asymmetry observed 
by easy axis hysteresis loop. In our samples the forward loop shows a slope, 
as if the full magnetization reversal takes place slowly, while the recoil loop is 
very square which may be the consequence of either coherent magnetization 
fiip (such as in monodomain particles) or the presence of high mobility domain 
walls. 



3 Model 



The magnetization reversal that occurs in exchange biased bilayers can not be 
explained with a simple algebraic addition of an exchange field to the applied 
magnetic field Ha because of the asymmetric shape of the hysteresis loop. 
In order to take into account the exchange bias, Meiklejohn and Bean [1] 
have proposed a linear term that accounts for the hysteresis curve shift, but 
which does not explain the asymmetry of the hysteresis loops for exchange 
biased F/AF systems. Owing to the large variety of hysteresis asymmetric 
curves, one should admit that one parameter exchange energy is no longer 
sufficient to interpret the form of these curves. In order to take into account 
the complex physical phenomena giving rise to the asymmetry of hysteresis 
curve, other terms should be added to the exchange energy expression describ- 
ing the F/AF exchanged biased samples. Several authors added a non linear 
biquadratic term proposed by Slonczewski[17]. However, that term is not suf- 
ficient to account for a large variety of asymmetric hysteresis curves known 
for various systems of F/AF samples. Based on symmetry considerations Am- 
brose et al. [6] proposed a Fourier odd and even series development of the shift 
field and the coercivity respectively that accounts for the angular variation of 
Hs and He- More recently, using the Stoner-Wohlfath model, involving only a 
uniaxial and a unidirectionnal exchange coupling, Xi et al. [10] fitted the hard 
axis hysteresis loop to experimental data. Nevertheless, the easy axis hystere- 
sis loop shape fitted poorly with experimental curve. Different mechanisms of 
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magnetization reversal such as, the AF domain wall formation, the different 
direction of easy axis in AF and F layers, the grain formation etc contribute 
to the hysteresis curve shape. Due to the complexity of magnetization reversal 
processes, which are strongly influenced by the physical state of F/AF interface 
and therefore very sensitive to the physical processes of sample preparation, 
a large number of Fourier series terms should be taken into account in order 
to achieve an agreement between the experimental data and theoretical cal- 
culations. In general, Fourier series convergence is rather slow. This is why we 
propose a cosine power series expansion of the exchange energy. Moreover, the 
grain formation and the easy axis distribution of different grains with respect 
to the magnetization direction of F layer leads us to introduce a mean easy 
axis direction of different AF grains. The exchange coupling between different 
grains, and the way they are influenced by F layer, suggests that statistical 
mean easy axis of AF grains may deviate from that of F layer. The exchange 
coupling with the AF layer is assumed to result in a single domain in the 
F film as recently observed in NiFe/Fe^oMn^Q by magneto-optic Kerr effect 
[18]. Therefore, we admit a uniform switching of the total magnetization of 
the F layer. The magnetic configuration of the F/AF bilayer is shown in Fig.l. 
The magnetic energy per unit area of exchange coupled F layer and AF layers 
can be written as follows: 

E = Kptpsin^Lp - HaMptFCOsie -if)- E^{(p - (3) (1) 



where Kp is the anisotropy constant of the F layer with its magnetization 
Mp making an angle ip with respect to the the F easy axis. The first term 
represents the anisotropy energy of the F layer. The second term corresponds 
to the Zeeman energy of the F layer submitted to an external static field 
applied at an angle 9 relative to the anisotropy easy axis of the F layer. The 
last term represents the exchange anisotropy energy, where j3 is the angle of 
the AF mean statistical easy axis with respect to the F easy axis. Ee{(p-P) 
may be written as: 

E,{<f-(3)^Y.JnCos^{f-P) (2) 

n 



If one admits, a priori, that Ec{f~f3) is linear: E(.{^p — (3)=JiCos{Lp — f3) then 
that represents the simple unidirectional energy originally proposed by Meikle- 
john and Bean [1] (with /? = 0) and used recently by Xi et al [10] (with /3 7^ 0). 
The biquadratic term of Ef.{(p — (3) was already used by Slonczewski [17] and 
later by Stamps [13] to describe the coupling of the F layer to both sublattices 
of two sublattice antiferromagnet. Sometimes, based on crystal symmetry con- 
siderations, threefold and fourfold components were used in Ee{'p ~ energy 
functional. The threefold component was recently found to be responsible for 
the symmetry breaking in the magnetization reversal process in Fe/MnF2[ll]. 
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Tang et al. [12] have shown that the exchange bias effect, described by a uni- 
directional anisotropy, is also accompanied by induced uniaxial and fourfold 
in-plane contribution in Fe/MnPd exchange biased bilayers. Because of the 
complexity of the exchange bias phenomenon, we can not limit our develop- 
ment to any finite number of terms. So the total energy in terms of effective 
field can be written as: 

E H 

— — = -^sin^cp - HaCos{e - (p) - V i7„cos''((^ - (3) 

IfMf I 



Where Hk is the anisotropy field and the are the coefficients of the power 
series development. Based on the above model, the magnetization curves of 
the F/AF bilayer are obtained by numerical calculations. 

The coefficients Hi, H2, H3 etc... being the power expansion coefficients of en- 
ergy functional, do not reflect the overall energy symmetry. We have truncated 
cosine power series at n= 5 for which the flt to the shape of easy and hard 
axis hysteresis curves was the best. The fitting to the easy axis hysteresis loop 
improves when one takes into account more terms. However, in that case, the 
hysteresis loop calculated for 9^{) by using fit parameters of easy axis hystere- 
sis loop shows parasitic minor loops. One may obtain better fit by using 10 or 
20 series terms. But that will not mean anyhow that the crystal has a 10 th 
or 20-th order symmetry since only the sum of all terms should provide such 
or such symmetry and not the each term. When we truncated the expansion 
to n = 4 the best fit to easy axis hysteresis curve becomes square. 

To give a physical signification to the 5 th order coefficient is difficult since 
our samples were prepared with RF sputtering and none of the textures were 
pronounced. After a thermal annealing, an interdiffusion have been induced 
[19]. Moreover, as reported in a previous work [15], a ternary NiMnFe alloy 
seems to be formed in the F/AF interface region. Chen et a/. [19] showed that 
the interdiffusion alters dramatically the exchange bias and that the hysteresis 
loop shape becomes asymmetrical. To the best of our knowledge there is no 
systematic study concerning on the infiuence of the alloy formed at the inter- 
difusion layer at the F/AF interface on the exchange energy functional. So it 
is a difficult task to associate a physical meaning with the exchange energy 
expansion coefficients. 



4 Results and discussions 



In a first step we have fitted the hysteresis loop to the theoretical model given 
by Eq. 3, when the field is applied along the easy axis. The comparison be- 
tween theoretical and experimental loops are presented on Fig. 2(a). One can 
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observe a strong asymmetry between the left and right side of the measured 
hysteresis loop. The magnetization reversal on the left side of the loop exhibits 
a hard-axis-type slope while the reversal for the right-hand side of the loop 
occurs with a vertical branch of the loop. It can be observed that the theo- 
retical magnetization curve fits well the measured hysteresis loop. The model 
gives values expected either for the hysteresis loop shift {Hs — 950 A/m) and 
coercivity (He = 900 A/m) or for the strong asymmetry in the magnetization 
reversal. The values of the fit parameters are /5=-0,064, Hk = 533 A/m, Hi 
= 1661 A/m, H2 =-903 A/m, =616 A/m, H4 =501 A/m, =-599 A/m. 
The first term is about two times larger than the anisotropy field of the 
as-deposited NiFe layer (i.e. 360 A/m). The second term Hi, that may be 
related to the exchange biased field, is about 1.5 larger than the hysteresis 
loop shift. 

One can observe a distinct kink on the left side of the calculated magnetization 
loop whereas the experimental data show a smooth variation. Leigthon et 
a/. [20] attributed that kind of kink to a two-stages magnetization reversal. The 
defect concentration is of high importance for exchange bias characterization 
as indicated by several authors (see the review by Stamps [13]). The higher 
order Fourrier coefficients in the exchange energy functional of single crystal 
samples are generally much larger than those for the poly crystalline samples. 
In polycrystalline materials the situation is complicated due to the spread 
in AF grain concentration. It is further complicated in the presence of an 
inter diffusion layer that may be a ternary NiMnFe alloy. According to Mewes 
et al. [21] the Stoncr-Wolfarth model shows sharp edges which arc rounded 
in the experiment for polycrystalline samples. This may explain the smooth 
variation of our experimental results. 

Energy diagrams are often useful to explain the magnetization reversal in 
magnetic hetorostructures. Fig. 3 represents the left (Fig. 3(a)) and the right 
(Fig. 3(b)) side of the measured hysteresis loop. In our calculations, for each 
value of the applied field, the system is in a local energy minimum which is 
plotted for the forward (left side) and reverse (right side) loops on the same 
figure. The energy curve presents an hysteretic behavior as a function of the 
applied field. On the same figure one can observe a strong asymmetry in the 
magnetization reversal due to the slow variation of energy local minimum as a 
function of applied field. On the right side of the figure (Fig. 3(b)), there is an 
abrupt lowering of local energy minimum at about Ha=0 which corresponds 
to an abrupt change of the energy minimum angle (ip). At the left side of the 
loop, the magnetization must overcome an energy barrier that vanishes after 
an increase AH of the field (Fig. 3(a)). This reversal may occur with a number 
of discontinuous jumps giving rise to "kinks" in the hysteresis loop. For the 
reverse loop, the magnetization vector fiips rapidly from (p=7i to ip=0 at a 
critical field about i?a— 0- 
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We have used the same fit parameters to calculate the magnetization curves for 
an applied field making different angles (vr/3, and 7r/2) with the F easy axis. 
The results are presented in Fig. 2(b and c). As observed on the Fig. 2( b and 
c) , the small deviation of experimental and theoretical curves is attributed to 
the initial dispersion of the as-deposited NiFe anisotropy revealed by a small 
coercivity mcasTircd along the hard axis. That may come from the aligning 
field in the sputtering chamber which lacks uniformity and may give rise to 
the F anisotropy dispersion. Indeed we did not include in our model such a 
dispersion term which is the main cause of that difference. Fig. 4 describes the 
comparison of theoretical and experimental angular variation of the coercivity 
He Fig. 4(a) and the hysteresis loop shifts Hg Fig. 4(b). The fit parameters 
correspond to the easy axis hysteresis loop. The behavior of the coercivity is 
very similar to that found by Ambrose et aL[6]. The measured coercivity does 
not vanish in the directions different from that of the easy axis (Oy^Q). For 
some angles there is discrepancies between the computed hysteresis shift and 
the measured one. First, the model predicts maxima of H,^ located at about 
6 = ±7r/15 and 7r±27r/3. This behavior is different from the believed one, 
according to which the largest values of Hg are expected at 6* = and tt as 
mentioned by several authors [6,10,14]. We observed experimentally the largest 
values of Hg ior 9 — and tt. Our angular dependence of the coercivity and 
hysteresis loop shift differs from the published results [6]. In our calculations 
the theoretical angular variation of He and Hg is obtained from a unique set 
of parameters deduced by the fit to the easy axis hysteresis loop instead of 
two different fits for the coercivity and the exchange field as usually presented 
[6,7]. 

We think that, the discrepancy between theoretical and experimental results 
in Fig. 4 is due to the fact that our samples arc polycrystallinc. Ricdiing et al 
[22] show that a complex 115(6*) behaviour may be due to the epitaxially grown 
samples. Liu et al [23] realised epitaxially grown NiFe/FeMn exchange biased 
bilayers on Si (100) films buffered by Cu. LEED experiments supported the 
existence of a six-fold symmetry. The obtained Il5(^) and IIc(^^) curves were 
strongly dependent on the growth technique. Polycrystallinc samples show 
perfect cos(6') behaviour even in the presence of a well pronounced induced 
(111) texture. That work supports our experimental data where also a well 
defined cos(^^) behaviour was observed. The epitaxial samples in the work of 
liu et al. show that H5'(^^) and Hc(^^) curves have several minima for different 
Cu buffer layer. The behaviour of l^s{d) and Hc(^^) curves shapes may also 
depend on the existence of an interdiffusion layer, but we have not sufficient 
experimental data to confirm or exclude such an influence. 

An important feature of the present model is the non zero value of /9=-0.064 
which is the angle between the F and AF easy axis directions. We suppose 
that it is due to the AF grains formation and the local exchange energy varia- 
tion at the interface of F/AF layers. Owing to the interface imperfections, the 
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easy axis of AF grains will suffer some dispersion [24,5]. A statistical mean 
easy axis direction can be a useful parameter to quantify that dispersion. The 
physics of such dispersion is somewhat complicated due to the lack of informa- 
tion about, the exchange energy local variation, the roughness characteristics 
and its influence on the energetic balance of concurrent energetic processes for 
energy minimization. This is why the value of j3, a priori, may be slightly dif- 
ferent from zero. So we think that it can be an additional pertinent parameter 
to describe the exchange bias. 



5 Conclusions 

In summary, we have presented a model based on the representation of the 
exchange energy by cosine power series. We have shown that it is possible 
to fit the asymmetric shape of the easy axis hysteresis loop. The asymmetric 
hysteresis loop shape is discussed on the basis of energy considerations. An 
angle f3 is introduced to quantify the easy axis dispersion of AF grains. The 
easy axis fit parameters describe satisfactorily the hysteresis curves in other 
directions. 
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Fig. 1. Vector diagram for an exchange coupled F/AF bilayer submitted to an 
applied field Ha 



Fig. 2. Representative hysteresis loops for a 

iVf8i-Fei9(520l)/Mn46A^i54(800A)biIaycr at (a) 0, (b) 7r/3, (c) it/2 of the 
applied field referred to the easy axis. The experimental data are shown by 
scattered empty dots. The solid lines are the theoretical magnetization curves. 



Fig. 3. (a): Forward measured magnetization hysteresis loop (upper) and corre- 
sponding calculated energy hysteresis loop (lower), (b): reverse measured magneti- 
zation hysteresis loop (upper) and corresponding calculated energy hysteresis loop 
(lower) . 



Fig. 4. Angular dependence of (a) shift field Hg and (b) coercivity of an exchange 
coupled NiFe/MnNi bilayer. The open dots are the experiments and the solid line 
is the theory. 
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Figure 4 

D. Spenato, S.P. Pogossian and H. Le Gall 

14 



References 

[I] W. Meiklejohn, C. Bean, Phys. Rev. 105 (1957) 904. 

[2] D. Mauri, H. Siegmann, P. Bagus, E. Kay, J. Appl. Phys 62 (1987) 3047. 

[3] N. Koon, Phys. Rev. Lett. 78 (1997) 4865. 

[4] A. Malozemoff, Phys. Rev. B 35 (1987) 3679. 

[5] M. Stiles, R. McMichael, Phys. Rev. B 59 (1999) 3722. 

[6] T. Ambrose, R. Sommer, C. Chien, Phys. Rev. B 56 (1997) 83. 

[7] N. Cheng, K. M. Krishnan, E. Girt, R. Farrow, R. Marks, A. Kellock, A. Young, 
C. Huan, J. AppL Phys. 87 (2000) 6647. 

[8] J.-V. Kim, R. L. Stamps, B. McGrath, R. Camley, Phys. Rev. B 61 (2000) 8888. 

[9] J. Geshev, L. Pereira, J. Schmidt, Phys. Rev. B 64 (2001) 184411-1. 

[10] H. Xi, R. White, J. Appl. Phys 86 (1999) 5169. 

[II] I. Krivorotov, C. Leighton, J. Nogues, I. K. SchuUer, E. D. Dahlberg, Phys. 
Rev. B 65 (2002) 100402-1. 

[12] Y. Tang, B. Roos, T. Mewes, A. Prank, M. Rickart, M. Bauer, S. Demokritov, 

B. Hillerbrands, X. Zhou, B. Liang, X. Chen., W. Zhan, Phys. Rev. B 62 (2000) 

8654. 

[13] R. Stamps, J. Phys. D:appl. Phys. 33 (2000) R247. 

[14] M. Stiles, R. McMichael, Phys. Rev. B 63 (2001) 064405 1. 

[15] D. Spenato, J. Ben Youssef, H. Le Gall, O. Ostorero, J. Appl. Phys. 89 (2001) 
6898. 

[16] J. Nogues, I. K. Schuller, J. Magn. Magn. Mat. 192 (1999) 203-232. 

[17] J. Slonczewski, Phys. Rev. Lett. 67 (1991) 3172. 

[18] S. Zhou, K. Liu, C. Chien, Phys. Rev. B 58 (1998) R14717. 

[19] M. Chen, C. Tsang, N. Gharsallah, IEEE Trans. Magn. 29 (1993) 4077. 

[20] C. Leighton, M.R.Fitzsimmons, P. Yashar, A. Hoffmann, J. Nogues, J. Dura, 

C. Majkrzak, L Schuller, Phys. Rev. Lett. 86 (2001) 4394. 

[21] T. Mewes, H.Nembach, M.Rickart, S.O. Demokritov, J.Fassbender, 
B.Hillebrands, hys. Rev. B 65 (2002) 224423. 

[22] S.Riedling, M. Bauer, C. Mathieu, B. Hillebrands, R. Jungblut, J. Kohlepp, 
A. Reinders, J. Appl. Phys. 85 (1999) 6648. 

[23] C. Liu, J. Du, J. Barnard, G. Mankey, J. Vac. Sci. Technol. A 19(4) (2001) 
1213. 

[24] E. Fulcomer, S. Charap, J. Appl. Phys 43 (1972) 4190. 



15 



